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Abstract 

We concider, the blow-up solutions for a coupled reaction diffusion system 
with gradient terms. The main purpose is to understand whether the gradient 
terms effect the blow-up properties. We derive the upper and lower blow-up rate 
estimates under certain assumptions. 

1 Introduction 

In this section, we consider the Cauchy (Dirichlet) parabolic problem: 

u t = Au+ \Vu\ qi +v Pl , v t = Av + \Vv\ q2 +u P2 in Ox(0,T), 
u(x, 0) = uo(x), v(x, 0) = vo(x), in Cl, 

where p±,p2,£ (1, oo), q\, q% € (1,2], uo,vo > are nonzero smooth bounded 
functions on ft (not necessarily radial), = R n or Br. Moreover, in case of 
Vt = Br, u,v are further required to satisfy the condition: 

u(x,t) = 0, v(x,t) = 0, on afix[0,T). (1.2) 

The problems of semilinear parabolic equations have been studied by many 
authors, for instance, consider the Cauchy (Dirichlet) problem of the semilinear 
heat equation: 

Ut = Au + u p , in Jlx(0,T), (1.3) 

where p > 1, CI = R n or Br. It is well known that every positive solution blows up 
in finite time, if the initial data is nonnegative and suitably large |6j, [9] . Moreover, 
it was proved in [5JIT7] that the blow-up rate estimate for (jl.3p takes the following 
form 

u(x,t) <c(T -t)~^, (x, t) G n X (0, T). 
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Later, in [10] it has been shown that if we add a positive gradient term to the 
equation (jl.3p . namely 

u t = Au+ \Vu\ q + u p , (1.4) 

then that enhancing blow-up, and the influence of the gradient term becoming 
more important as the value of p decreases. In the case q = 2 for radial positive 
solutions in R n , it was shown in [8] that blow-up solutions behave asymp- 
totically like the self-similar solution of the Hamilton- J acobi equation without 
diffusion (ut = \S/u\ 2 + u p ), which takes the form 

u(x, t) = (T- t)^w( * ), m = (2 - p)/2(p - 1), 

where w 6 C 2 (R n ) is a positive radial decreasing function. On the other hand, 
the existence of nonnegative global solutions is shown in [15] for small initial data. 

In [3l [3j, it was considered, the Cauchy (Dirichlet) problem of the following 
semilinear system: 

u t = Au + v p \ v t = Av + u P2 , (x,t) E n x (0,T), (1.5) 

where p±,P2 > 1, ^ = Br or i? n , with nonzero initial data uq,vq > 0, it was 
shown that any positive solution of this problem blows up in finite time if the 
initial data are large enough. Moreover, for the Cauchy problem of (jl.5p . it is 
well known [3] that any nontrival positive solution blows up in finite time, if 

max{a,/3}>-, (1.6) 

where 

a = -, = -. (1.7) 

pxp 2 - 1 P1P2 - 1 

The blow-up rate estimates of this system was studied in [TJ [2], it was proved 
that there exist a positive constant C such that 

u(x,t) < C{T-t)~ a , (x,t) eCl x (0,T), 

v(x,t) < C(T-t)- p , (x,t) G CI x (0,T). 

In this paper, for problem (jl.ip . under some restricted assumptions, we prove 
that the upper blow-up rate estimates of the positive solutions and their gradient 
terms, take the following forms: 

2(P1+1) 

u(x,t) + |Vu(x,t)|MP2+2pi+i <Ci{T-t)~ a } (x,t) S x (0,T), 

2(P 2 + 1) „ 

u(x,t) + |Vu(x,t)|pif2+2p 2 +i <C 2 {T-t)~ li , (x,t) e x (0,T), 
where Ci, C 2 > 0. 
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2 Preliminaries 



Set 

Fx(v,Vu) = \Vu\ qi +v pi , F 2 (u,Vv) = \Vv\ q2 + u P2 . 

Since the system (jl.ip is uniformly parabolic and Fi,F 2 are C 1 ([0, oo) x R n ), 
moreover, the growth of the nonlinearities F± and F\ with respect to the gradient 
is sub-quadratic, it follows that, the local existence of the unique nonnegative 
classical solutions to the Dirichlet (Cauchy) problem of (jl.ip is guaranteed, for 
smooth and bounded initial data uq,vq,, by the standard parabolic theory [11] 
(see also |12j). On the other hand, the positive solutions of problem (11. ip may 
blow up in finite time, and that due to the known blow-up results of the system 
(I1.5P and the maximum principle |13j . 

Remark 2.1. Since the growth of the nonlinear terms in problem (jl.ip with 
respect to the gradients is sub-quadratic, the gradient functions Vn, Vv are 
bounded as long as the solution (u, v) is bounded (see [12]). 



2.1 Blow-up Rate Estimates 

In the next theorem, we establish the upper blow-up rate estimates for the prob- 
lem (II. ip . Furthermore, without comparing the blow-up solutions of this problem 
with those of problem fjl ,5j) , we show that the blow-up can only occur simultane- 
ously. 

Theorem 2.2. If Pi,p 2 ,qi and q 2 satisfy the following conditions 

(1) max{a,/3} > §, 

(2) k i<s 2 <!S, 

where a, f3 are given in |-?.7| ), then for any positive blow-up solution (u, v) of 
problem there exist positive constants Ci,C 2 such that 

2(P1 + 1) 



u(x,t) + |Vn(x,t)|pif2+2Pi+i <d(T-t)~ a , (2.1) 

v(x,t) + |Vi)(x,t)|fif2+2P2+i <C 2 (T-t)~ (j , (2.2) 
in Q x (0, T), where T < oo is the blow-up time. 
Proof. Let 

M u (t) = sup [u{x,t) + \Vu{x,t)\piP2+ 2 Pi+ 1 ], 
nx(o,t] 

2(P2+1) 

M v (t) = sup [v(x,t) + |Vv(x,i)|pira+2p2+i], 

f7x(0,t] 
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for t G (0,T). 

Clearly, M u , M v are positive, continuous and nondecreasing functions on (0, T). 
At least one of them diverges as t — > T, due to (it, t>) blows up at time T. 
We show later that there is 5 G (0, 1) such that 

5 < M~^(t)M^(t) < i t g (T/2,T). (2.3) 

o 

So that, consequently, both M U ,M V have to diverge as t — >■ T. 

To establish the blow-up rate estimates, we use a scaling argument similar as 
in pp. The proof is divided into several steps. 
Step 1: Scaling 

If M u diverges as t — > T, the following procedure can be applied. 

Given £ G (0,T), choose (x*,t*) e(]x (0,i ] such that 

2(P1+1) I 

u(s*,t*) + \Vu(x*,t*)\piP2+2 P1 +i > -M u (t ). (2.4) 

i 

Let 7 = 7(£o) = M u 2a (to) be a scaling factor. Define the rescaled functions 

<pl(v,a) = l 2a u{ 1 y + x*, 1 2 s + t*), (2.5) 

<4{y, s) = i^v(ciy + x*, 7 2 s + **), (2.6) 
for (y,s) E0 7 x (-7~V, 7~ 2 (T - t*)), where 

fi 7 = {y G R n : jy + x* G fi}. 

Clearly, 



# n if n = # n , 

Bii if n = B R . 



Next, we aim to show that (tpj, tp^) is a solution of the following system 



where /zi = 2a + 2 - (2a + /i 2 = 2/3 + 2 - (2/3 + l)g 2 . 
From the assumption (2), it follows that > 0. 

Clearly, 

^= 7 2 «+ 2 u, V^ = 7 2o+1 Vn, A^= 7 2a+2 An. (2.J 
From P^]> - it follows 



Multiply the last equation by ^ 2a+2 \ we get the first equation of (|2.7p . In the 
same way we can show that ^ satisfies the second equation in system ()2.7[) . 
Restrict s to s 6 (— 7 _2 i*, 0], our aim now is to show that 

2(pi+1) 

?7(y,s) + |V^(y,s)|^2+2P 1+ i < l, (2.9) 

for (y,s) e Oy x (-7~ 2 t*,0]. 
Prom (12,81). we obtain 



2(P1+1) , 2( P1 + 1) 2( P1 +1) , 2( P1 +1) 

[Vy7(y, S)| P1P2+2P1 + 1 = -y L PlP2- 1 JI 'PlP2+ 2 Pl+ lJ |Vn| p l p 2+ 2 Pl+ 1 , 

2( P1 + 1) 

= 7 2a |V'u| J, i p 2+ 2 Pi+ 1 . (2.10) 



Clearly, 



2(P1+1) 

«(x,t) + |Vu(a;,i)l P1P2+2pi+1 < M tt (t ), 0, t) € O x (0, t*}. (2.11) 

From ([23]), (|2~10D and flgZEED , we get (1231) . 

Moreover, 

2( P2 + 1) _£ 

^ + |V^|^2+2p 2 +i < M u » (t )M v (t ), (2.12) 

for (y,s) £li 7 x (-7~ 2 t*,0]. 

On the other hand, from (]2.4h . we obtain 

2( Pl +l) i 

¥>I(0,0) + |V^(0,0)|fif2+2pi+i > - (2.13) 

If M v diverges as t — > T we can proceed in the same way by changing the role of 
u and v. 

Step 2: Schauder's estimates 

We need interior Schauder's estimates of the functions ipi,<p2 on the sets 

S K = {y 6 n 7 , Is/I < if} x [-K,KL], K > 0, L = 0,1. 
Assume that </?i , satisfy in the condition 

2( P i+l) 2( P2 +1) 
< (f[ + |V^|P1"2+2P1+1 < £?, < + |V^|W2+2P2 + l < 5. (2.14) 

We claim that for any K > 0, B > and <r > small enough, there is a constant 
C = C{K,B,a) such that 

lKH c2+CT , 1+ § (Sjf) < C, ||^|| c 2+.,i +f (Sjf) < C. (2.15) 

From (|2.14p we deduce that each of V?i > V2 > V^pJ , V<^2 > is uniformly bounded 
function in S2K- Therefore, the functions ((pj) pi , ((p1) P2 , |V(^| 91 , (V^l 92 are uni- 
formly bounded in S2K- So, the right hand sides of the two equations in (|2.7p are 
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uniformly bounded functions in S2K, applying the interior reqularity theory (see 
we obtain (locally) uniform estimates in C 1+<T,_ 2 _ -norms. Consequently, 
by Lemma ??, we obtian (locally) uniform estimates in Holder norms C 7 '? on 
the right hand sides of the both equations in (|2.7p .Thus the parabolic interior 
Schauder's estimates imply (|2.15p (see [?. [TT]). 
Step 3: The proof of ([23]) 

Suppose that this lower bound were false. Then there exist a sequence tj — > T 
such that 

M u 2a (tj)M v 2p (tj) — ► 0, as j -> 00. (2.16) 

Then clearly M u diverges as tj — > T. For each tj in the role of to from Step 1, we 
scale about the corresponding point (xptj) for all j such that t* < tj, with the 
scaling factor 

lj = 1 (tj) = M~^{t 3 ). 
We obtain the corresponding rescaled solution (cfi , <p!i)i 

¥>? (V, s) = ^ a u( ljy + x*, 7 ]s + t*), (2.17) 

$ (V, ») = 7f vfrjV + x*,^s + t*j). (2.18) 
Clearly, (ip^ , <p^ ) satisfies (as in Step 1) the following problem 



^-A^= 7 f|V^^ + (^r, 

with 



for (y,s) G 7j . x (-7^. 2 t*,0], where 




if fi = J R n , 
if Q = B R . 



Clearly, 

— > R n , as j — > 00. 

From (I2~TH1) . flZZQj) , we see that 

2 (P2 + !) 

+ | P1P2+2p2+1 — > 0, as j — )■ 00. 

Thus ,V<P2 are bounded in f] 7j x (— 7~ 2 i*,0], Vj. 



(2.19) 



2( P1 +1) 

pj j (0,0) + |V^ j (0,0)|w2+2pi+i > 1/2, 

< Pi + |V<^'|piP2+ii+i < 1, ^ (2-20) 

2(p 2 +l) _£ 
^ + I I P1P2+2P2 + 1 < M M a (tj)M v (tj), 

2.* 
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Using the uniform Schauder's estimate derived in Step 2 to {ip^ , ip 2 3 ) 
H^H^+a.i+f ({ y eQ^.,\y\<K}x[-K,o]) - ° K > 

||^ J || c2+CT , 1+ f ( {ye n 7 .,|j,|<^}x[-^,o]) - Ck > 

where Ck is independent of j. 

Since (cpj 3 ', ip 2 3 ) is defined on a compact set, by the Arzela-Ascoli theorem, 
there exist a convergent subsequance, we still denote it by (ipj 3 , tp], 3 ) ■ 

Since fii,fJ,2 > and V^, VyJ are bounded, it follows that, the limit point 
(yii ^2) is a solution of the following system 

cp ls = Aw + tpf , ip 2s = A^ 2 + <Pi\ in J R n x(-oo,0]. (2.21) 

Since (p 2 3 — > 0, as j — > 00, it follows that 992 = 0, in R n x (— 00, 0]. 
Consequently, from the second equation in (|2.2ip . we obtain that 

991=0, in i? n x(-oo,0]. 

This means 

2(P1+1) 

99i(0,0) + |V99i(0,0)|wra+ 2 Pi+ 1 = 0, 

which contradicts with (|2,20p . Thus, the lower bound is held. 

To prove the upper bound of (|2.3|) we proceed similarly as in the proof of 
lower bound with changing the role of u and v. 
Step 4: Estimate on doubling of M u 

As M u is continuous and diverges as t — > T, for any to € (0, T) we define t$ 

by 

t+ = max{t e (t ,T) : M u (t) = 2M u (t )}. 

Clearly, 

2(P1+1) 

u(x,t) + \Vu(x,t)\piP2+2Pi+i < 2M u (t ), (x,t) e!lx {0,t+}. (2.22) 

Take 7 = 7 (to) = M~^(t ). 
We claim that 

7~ 2 (^o)(*o ~~ *o) < A to€(-,T), 

where the constant A £ (0,oo) is independent of to. Suppose that this estimate 
were false, then there would exist a sequence tj — > T such that 

l~\tj)(t+ -tj)^ 00, 

where 

tt = max{t £ (^ , T) : M u (t) = 2M u (tj)}. (2.23) 
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For each tj we scale about the corresponding point (Xj,t*) such that 

0<t*<tj, ^ <tj <t+<T, Mj 
with the scaling factor 

lj = 1 (t j ) = M~^(t J ). 

As in Step 3, we obtain the corresponding rescaled functions (ip^ , ip^ 3 )> which 
satisfies (|2.19p with the following conditions 



2( P1 + 1) 

</>i J '(0,0) + |V^ J '(0,0)|fif2+2pi+i > 1/2, 

< + |V^ |P1P2 + 2P1 + 1 < 2, 

2(p 2 + l) _§_ 
$ + \Vifl j I P1P2+2P2 + 1 < M u <* (tj)M v (tf), 



(2.24) 



for (y,s) G 7 . x (-7. V, 7 . 

From (|2~23jl and (l234j) . it follows that 



2(p 2 +i) a _ £ 

^ + |V<^>iP2+w < 2«M„ Q (tt)M„(tt). (2.25) 



From (12.31). we have 



M v (t) < 5- 2 ?M«(t), te(^,T). 
Therefore, (|2.25[) becomes 

y , -■ , 2 <P2 + 1) 2« 

^ + |V<^ | P1P2+2P2+2 < • 

By using the Schauder's estimates derived in Step 2 for (ipj 3 , iplf ) we get a con- 
vergent subseguence in C^ c a ' l+<T ^ \R n x R) to the solution of system (|2.2ip in 
R n x R. This is a contradiction because all the nontrival positive solutions of 
system <|2.21 j> . under the assumption (1), blow up in finite time (see [3]). 
Thus, there is A > such that 

T 2 (to)(4-t ) < A, t oE (|,T). (2.26) 

Step 5: Rate estimates 

As in Step 4, for any t G (T/2, T) we define 

t x = £+ G (t ,T) such that M u (fi) = 2M„(t ). 
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Due to (gJSD , 

(h -t ) < AM~"{t ). 
We can use t\ as a new to and obtain £2 G (£> T) such that 

M„(t 2 ) = 2M u {t 1 ) = 4M u (t ), 

(*2 -ti) < Afcf« "(ti) = 2-aAM^ s («o). 
Continuing this process we obtain a sequence tj — > T such that 

fo+i " *j) < 2"- AMJ" - (to), j = 0,1,2,... 
If we add these inequalities we get 

(T-t ) < J2 2 ~" AM u*(to)- 

j>0 

Thus 

(T - t ) < (1 - 2-^y 1 AM~« (t ) 
From using (|2.3p we obtain 

M w (to) < 5~ 2t3 M§(to), t G (T/2,T). 

Thus 

M v (t ) < 5- 2 P(l - 2-i)-PA?(T - to)~P, t G (T/2,T). 
From above there exist two constants C* , C 2 such that 

M u (to)<C\(T -toY a , t e(^,T), 

M v (to)<C* 2 (T-t )-P, *oG(|,T). 

From the last two equations and the definitions of M U ,M V , it follows that there 
exist constants C\ , C 2 such that 

2(P1 + 1) 

u(x,t) + |Vn(x,t)|pif2+2Pi+i < Ci(T-i)" Q , 

2(p 2 + l) 

u(cc,i) + |Vv(x,t)|piP2+2p 2 +i < C 2 (T - t)~P, 
for (x,t) G X (0,T). □ 
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Remark 2.3. If no = i>o, P = Pi = P2-, Q = Qi = 12, then problem (jl.ip can be 
reduced to a scalar Dirichlet (Cauchy) problem for (jl.4p . Moreover, if 

2 2r> 
l<p<l + _, i< g <_^, (2.27) 
n 1 + p 

then in a similar way to the proof of Theorem 12.21 we can show that, for a 
nontrivial positive blow-up solution u, there exist C > such that 

u(x,t) + \Vu(x,t)\& < C(T-t)t=*, in nx(0,T), (2.28) 

i.e. 

u(o:,t) < C(T-t)i=r, in Ox(0,T). (2.29) 

As we have mentioned before, the rate estimate (|2.29|) is also known for the 
blow-up solutions of equations (j 1 .3 j) . Therefore, if p, q satisfy (|2.27p . then the 
positive gradient terms which appear in equation (|1.4|) . does not affect the blow- 
up rate estimates of these problems. A similar observation holds for problem 
fll.l|> by Theorem 12.21 which shows that the upper rate estimates of the Cauchy 
or Dirichlet problem for system (jl.ip are the same as those known for the system 
(11. 5p . Therefore, under the assumptions of Theorem 12.21 the gradient terms in 
system (jl.ip have no effect on the blow-up rate estimates. 

2.2 Blow-up Set 

It is well known that for the semilinear system (|1.5p defined in a ball, under some 
restricted assumptions on uq, vq (nonnegative, radially decreasing functions), that 
the only blow-up point is the centre of that ball (see [33]), while it is unknown 
whether this holds for the system (jl.ip . However, for the radial solutions of 
the single equation (|1.4j) defined in fi, in case q = 2, there is global blow-up, if 
1 < p < 2, £2 = Br or R n , and regional blow-up, if p = 2, $7 = R n , while a single 
blow-up point, if p > 2, $7 = Br (see [131 IS] an d the references therein). The 
proof relies on the transformation v = e u — 1, which converts fll.4|> into the 
semilinear heat equation vt = Av + (l + v) log p (l + f). We note that, these results 
are much different from those known for equation (jl.3j) (see [H]), because for 
any p > 1, fi = Br or R n , only a single blow-up point is known to occur for 
that problem, where the initial date are nonnegative, radially nonincreasing and 
bounded function. 
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